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Abstract. We consider some recent developments in GRB/afterglow observations: (i) the appear- 
ed ■ ance of a very hard prompt component in GRB 941017, and (ii) variability in the early afterglow 
q | light curves of GRB 021004 and GRB 030329. We show that these observations fit nicely within 

the internal-external shocks model. The observed variability indicates that the activity of the inner 
^"^ ■ engine is more complicated than was thought earlier and that it involves patchy shells and refreshed 

shocks. We refute the claims of Berger et al. and of Sheth et al. that the radio and mm observations 

of GRB 030329 are inconsistent with refreshed shocks. 

^ ■ Introduction 

The early afterglow and the GRB/afterglow transition are among the unexplored 
regimes of GRBs. A great progress on this front was achieved during the last year 
when, following fast HETE II identifications, two afterglows (GRBs 021004, 030329) 
were followed from very early on showing remarkable variability and rich structure. 
Also, somewhat unexpectedly, a search within the BATSE/EGRET archives revealed a 
new very hard long lasting (~ 200 s) component of GRB 941017. This is most likely 
a manifestation of the early afterglow and of the GRB/afterglow transition. We discuss 
these developments and their implication within the internal-external shocks model. 
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The Prompt High Energy Emission from GRB 941017 as a 



GRB/afterglow transition 

Recently, Gonzalez et al. |Q]] discovered a high energy tail that extended up to 
200 MeV in the combined BATSE and EGRET data of GRB 941017. The tail had a hard 
spectral slope (F v oc v°) up to 200 MeV. It appeared ~ 10-20 s after the beginning of the 
burst and displayed a roughly constant flux while the lower energy component decayed. 
At late times (~ 150 s after the trigger) the very high energy (~ 10-200 MeV) tail had a 
luminosity ~ 50 times higher than the "main" y-ray energy band (~ 30 keV-2 MeV). 

Granot & Guetta [2] where the first to suggest that we see here another manifestation 
of the very early afterglow. Sari |3J has shown that for long bursts the external shocks 
(and hence the afterglow) begin ~ R/cT 2 after the beginning of the burst while the 
internal shocks are still going on and hence the burst is still active. This behavior 
was^seen in the transition from the harder initial burst to the softer early afterglow 
It was also seen in the prompt optical flash of GRB 990123 where the lower 



(optical) energy component did not trace the ~MeV y-rays and a pronounced hard to 
soft evolution was seen in the y-ray signal. 

The very high energy of this emission suggests that it is inverse Compton. There 
are two relevant emitting regions in the early afterglow: the forward shock and the 
reverse shock. With typical parameters $7$ the energy of the synchrotron photons and 
the electrons' Lorentz factor in the forward and reverse shocks are: 

v sy nch,F « 0. 1 MeV (r/300) 4 ; y e . F « 10 4 (r/300) 
v sy nch,R « 1 eV (r/300) 2 ; y e , R w 300 

There are four possible combinations of seed photons and scattering electrons: 

Rev. shock electrons For. shock electrons 

Rev. shock photons ~ 0. 1 MeV(r/300) 2 ~ 100MeV(r/300) 3 

For. shock photons ~ 10GeV(r/300) 4 ~ 10TeV(r/300) 5 (within Klein Nishina) 

While these approximate results are very sensitive to T, they indicate that inverse Comp- 
ton scattering of the reverse shock photons on the forward shock electrons yeilds the 
right energy range. The detailed calculations of Pe'er & Waxman confirm these naive 
estimates. The main problem, however, is not to explain the location of the spectral peak 
but to explain the spectral slope (F v V ) and the temporal slope (F v °= t ). Pe'er & Wax- 
man [8] reproduce the spectral slope by requiring that the synchrotron self absorption 
frequency of the reverse shock emission would be high enough to effect the observed 
spectrum. Granot & Guetta B2I] reproduce both the temporal and spectral behaviors by 
considering a slightly different scenario where the high energy component is produced 
by Synchrtron self Compton within the reverse shock. They require a slightly higher ex- 
ternal density with a somewhat unusual profile, R~ l , for a uniform ejecta shell (which 
may explain the rareness of the event). In both models the high energy component is a 
clear manifestation of the onset of the afterglow and the GRB/afterglow transition. 



Afterglow Light curve variability 

Theory^ The different scenarios that lead to afterglow temporal variability can be 
distinguished according to their characteristic features. Density variations produce only 
weak fluctuations above the cooling frequency, V c , and cannot produce sharp changes 
in the light curve. Energy variations produce variability both above and below v c , and 
can arise either due to refreshed shocks or due to a patchy shell structure. These two 
mechanisms produce very different light curves. While the former produce a step-wise 
increase in the light curve, the later produces random fluctuations with a decreasing 
amplitude. 

Rees & Meszaros 0, Kumar & Piran ifioll and Sari & Meszaros ifllll suggested 
that slow shells take over the slowing down matter behind the afterglow shock and 
produce refreshed shocks. Slow shells with T s emitted from the source right after the 
fast ejecta, catch up and collide with the slowing down ejecta at an observer time 
t ~ 0.25(r. s /10) _8//3 (£'i S o,52/«o) 1//3 days (where E- lso ^2 is the isotropic kinetic energy 



in units of 10 52 ergs and hq is the external density in cm~ 3 ) when the ejecta's Lorentz 
factor drops slightly below Y s . The clearest feature of refreshed shocks is a monotonous 
increase in the overall energy. Therefore the observed flux can only increase (relative 
to the expected decay). The light curve has a step wise form with each step produced 
by the arrival of a single shell. This step wise structure is seen both above and below 
the cooling frequency with a similar amplitude. Each step (in the optical light curve) 
should be accompanied by a flare in low frequencies that is produced by the reverse 
shock which propagates back into the slow shell iflflLfllll . The time scale, At of the steps 
and the corresponding flares depends on their timing relative to the jet break. Before 
the jet break the refreshed shocks are "locally" spherically symmetric and therefore the 
angular time imposes At ~ t II 1 CXI . The intensity of the reverse shock flare in this regime 
is calculated in H 1 Ofl , and the decay after the peak is t~ 2 IT]]. In the post-break case 
the cold slow shells may not expand sideways (if cold enough). Then they keep their 
original angular size, (?,, which is smaller than r~ . In this case At ~ tj < t, where tj is 
the jet break time [12] and the transition is fast. A reverse shock flare is expected in this 
case as well. However, its frequency, intensity and temporal decay (which is expected to 
be steeper than in the spherical case) are much harder to calculate. 

Kumar & Piran II 1 3M suggested, in the Patchy shell model, that the shells have an 
intrinsic angular structure. As the blast wave decelerates the angular size of the observed 
region (~ 1/T) increases. The effective (average) energy of the observed region and 
hence the observed flux, relative to the expected decay, varies with time depending on 
the angular structure. The variability time scale is At ~ t II14I1 . The averaging over a 
larger and larger random structure leads to a decay of the envelope as ?~ 3 / 8 fl4l[la1 . An 
important feature of this scenario is the break of the axial- symmetry and therefore the 
production of a linear polarization. The variation of the polarization, both in degree and 
in angle, are correlated with the light curve variations Ii4,[3]. The variability will be 
observed both above and below the cooling frequency v c with a similar amplitude. 



Wang & Loeb 1 17], Lazzati et al. [18] and Nakar et al. lfl5ll considered External 
density variations. Such variations may result from ISM turbulence or from a variable 
pre-burst stellar wind. Wang & Loeb [17] analyzed the light curve resulting from mild 
density fluctuations due to ISM turbulence. They show that these density fluctuations 
can produce short time scale (At < O.lt) and low amplitude (~ 10%) fluctuations in the 
light curve. Lazzati et al. lflal . Nakar et al. lfl5ll and Nakar & Piran lfl9ll considered 
large amplitude spherical density fluctuations. A basic feature of the resulting light 
curve that distinguishes it from energy variations is that in the former the light curve 
is different above and below the cooling frequency, V c . Density variations produce only 
weak fluctuations above V c . The amplitude of the fluctuations above v c is at most tens 
of percents and it is much smaller than the amplitude of the fluctuations below V c lfl9ll . 

A second feature of density fluctuations is their inability to produce a sharp variation 
(either increase or decrease) in the light curve lEoTl . First, we note that because of angular 
spreading, spherical density drops cannot produce decays sharper than f ~ 2,6 lll9l.l2lll and 
even this decay is reached very slowly. More interesting is the fact that even a sharp 
density enhancement cannot produce a steep increase in the light curve. The earlier 
calculations lfl5i[l8l.[l9ll assumes that the ejecta can be described by a Blandford-McKee 
solution whose density profile varies instantaneously according to the external density. 



These calculations do not account, however, for the reverse shock resulting from density 
enhancement and its effect on the blast-wave. Thus the above models are limited to 
slowly varying and low contrast density profiles. Now, the observed flux depends on 
the external density, n, roughly as n 1 / 2 . Thus, a large contrast is needed to produce a 
significant re-brightening. Such a large contrast will, however, produce a strong reverse 
shock which will sharply decrease the Lorentz factor of the emitting matter behind 
the shock, T s h, causing a sharp drop in the emission below v c and a long delay in the 
arrival time of the emitted photons (the observer time is r , ). Both factors combine 
to suppresses the flux and to set a strong limit on the steepness of the re-brightening 
events caused by density variations. Note that while non spherical density fluctuations 
may lead to a steeper decline they usually do not lead to a steeper increase in the flux. 

Implications: The early afterglow of GRB 021004 showed clear deviations from 
a smooth power law decay, lasting from 0.04 days to 3 days. The fluctuations in the 
light curve where accompanied by fluctuations both in the degree and in the angle of the 
polarization (LI 4] and references therein). 

The steep decays after each bump imply that the variations do not result from re- 
freshed shocks. Thus, variable external density variations lfl5l HrL I22I1 and the patchy 
shell model were considered as possible explanations. Unfortunately, the X-ray 

observations are not detailed enough to clearly distinguish between density and energy 
variations (although the former are favored by W22I1 ). However, the sharp decays cannot 
be produced by "locally" spherical density variations [19] . Furthermore, the first bump 
requires, using the instantaneous Blandford-McKee approximation, an increase in the 
external density by a factor of ~ 10 over AR/R ps 0.05 lll8L ll9ll. Such a density contrast 
produces a mildly relativistic reverse shock which reduce by a factor of « 2, making 
the approximation inconsistent. Preliminary results lEoll of detailed numerical simula- 
tions (including both hydrodynamics and synchrotron radiation) suggest that this bump 
cannot be produced by density enhancement (due to the suppressed forward shock emis- 
sion, caused by the reverse shock). These results leave the patchy shell as the only viable 
explanation. Indeed, Nakar & Oren IU4I1 show that patchy shell can reproduce the light 
curve (including the sharp rise and steep decay of the first bump). They show that angu- 
lar energy profiles which produce the observed light curve produce also a polarization 
curve that fits the observed polarization. We conclude that angular energy fluctuations 
are the dominant process that produce the observed fluctuations in GRB 021004. 

In addition to the remarkable supernova signature, the optical afterglow GRB 030329 
has shown also a unique variability (112D and references therein). Several step- wise 
bumps, at t = 1.5, 2.6, 3.3 and 5.3 days, are seen after the jet break at tj ~ 0.5 days. The 
first bump was the largest and best monitored, but even with the less dense monitoring of 
the later bumps the step-wise profile (where after each bump the original decay slope is 
resumed) is clear. All the bumps had a short rise time At f» 0.4-0.8 d < t. The step-wise 
profile seems like a clear signature of post-break refreshed shocks 1 12], where At mtj <t 
because the later slower shells did not expand sideways before colliding with the faster 
earlier ejecta. Moreover, the energy injected in these shocks is 10 times the energy in 
the original blast- wave. This late (or rather slow) energy injection explains an additional 
peculiarity of this GRB: the low energy output in y-rays and in the early X-ray afterglow. 

Berger et al. 112311 suggested that the first bump and the energy deficiency can be 



explained by a two component jet: A slow and energetic component with a wide half- 
opening angle (17°) dominates the afterglow after 1.5 days, and a fast component with 
a narrow half-opening angle (5°) dominates the afterglow before 1.5 days. The slow 
component is observed only after 1.5 days^ ? dec « 0.5(r s /10) _8 / 3 (£'i SOi 52/no) 1//3 days, 
since only at this time its reverse shock consumes the slow shell. However, this model, 
which received a great publicity, predicts that the rise time, At, of the first bump should 
be of the same order as the observed time, f^ec- Furthermore, it predicts a smooth light 
curve after the first bump. Both predictions are contradicted by the observations. 

Millimeter observations of GRB 030329 JH show that at 100 & 250 GHz the flux is 
rather constant during the first week. Most surprising are two measurements at 100 GHz, 
one before the first bump (0.6-1 d) and one after (1.7-1.9 d), which show a constant 
flux. Sheth et al. |24] and Berger et al. [23] claim that these results support the two 
component jet model and reject the refreshed shocks model due to the lack of radio 
flares. However this analysis overlooks the fact that the encounter of the slow component 
in the two-components jet model with the external matter produces a reverse shock 
which should produce a radio flash. This flash is analogous to the optical flash produced 
by the deceleration of the fast component ED, and to the radio flare expected in the 
refreshed shocks model. The timing of this flash is /dec = 1.5 d and its magnitude is 
easily calculated. The contribution from this reverse shock at f^ec = 1-5 d is larger by a 

5/3 

factor of up to T s ' ~ 15 than the flux from the forward shock. Thus a very bright and 
fast fading millimeter flash is expected in the two component jet model. A more detailed 
calculation (using 111 ill and the parameters of the wide jet presented in H23lT ) shows that 
both V a and v m of the reverse shock at t^ ec are around 100-200 GHz and that the flux at 
V m is ~ 500 mJy which is an order of magnitude larger than the expected flux from the 
forward shock and than the observed fluxes at 100 & 250 GHz (~ 50 mJy) at this time. 

The main argument of lE3ll in favor of a two-components jet is the existence of a sec- 
ond jet break in the radio after tji ~ 10 d. However, even this argument is not strongly 
supported by the data. According to this model the flux below V OT should rise as t 1 / 2 
before tj2 and decay as z -1 / 3 after tj2- However, the data of W24I1 contradict this pre- 
diction. At 100 & 250 GHz the flux is constant before the passage of V m and the decay 
after this passage (at t = 6 & 8 d < tj 2, respectively) is steeper than t~ in in both bands. 
This looks like a clear signature of a post break behavior in the radio at t > 1 d. Now 
the radio observations at lower frequencies (< 22 GHz) do not conform with the simple 
post-break model (F v oc ? -1 / 3 ). The flux at these frequencies rises with time before the 
passage of V m . Interestingly enough, this radio behavior is exactly the one predicted as 
the post break radio behavior by [25], using a 2D relativistic hydrodynamical simula- 
tions. We find the striking similarity between the radio observations and Fig. 2 of Il25ll 
as a very strong support that the broad band data are totally consistent with a single jet. 

Sheth et al., Il24ll emphasize the fact that radio flare H1QI1 was not detected during 
the first bump at t ~ 1.5 days and argue that this rules out the refreshed shocks model. 
However, in the post-break refreshed shocks model, the radio flash is expected to be 
fainter and to decay faster than the radio flash in the two-components jet model (due to 
the lower energy and the lateral spreading of the slow shell as opposed to the more 
energetic and "locally" spherical wide jet). Thus, by assuming spherical symmetry, 
Sheth et al., [24] over estimate the expected flash in the post-break refreshed shocks 



scenario. It may be that the flash was missed by the sparse measurements due to its 
lower intensity and faster decay. A detailed (and highly non trivial) calculations should 
be done in order find out. 



Summary 

With an increasing flow of new observations we discover that GRBs and afterglows 
are richer than what was previously thought. The simple spherical theory had to be 
modified, first with jets and now with additional angular structure (patchy shells) and 
more extended velocity structure (refreshed shocks). The simple synchrotron theory has 
to be modified with inverse Compton scattering. One can worry, are we adding epicycles 
trying to revive a wrong theory? We don't believe so. Complications and variation 
are common in astrophysics and are found everywhere in nature. Moreover, the three 
main themes that have been introduced here: patchy shells, refreshed shocks and inverse 
Compton, were not invoked aposteriori to explain the new observations. On the contrary, 
all three have been suggested long ago. It is just natural and even reassuring to discover 
them when better data become available. 

The research was supported by US-Israel BSF. 
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